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■ ABSTRACT 
l> '■ 

' Context. Sub-millimeter and Far-IR observations have shown the presence of a significant amount of warm (few hundred K) and 

' dense (n(H2) > lO'* cm"'') gas in sources ranging from active star forming regions to the vicinity of the Galactic center. Since the 

1 1 1 main cooling lines of the gas phase are important tracers of the interstellar medium in Galactic and extragalactic sources, proper and 
, detailed understanding of their emission, and the ambient conditions of the emitting gas, is necessary for a robust interpretation of the 

a observations. 
Aims. With high resolution (7" - 9") maps (~ 3 x 3 pc-) of mid-7 molecular lines we aim to probe the physical conditions and spatial 

(~| I distribution of the warm (50 to few hundred K) and dense gas («(H2) > lO"" cm"') across the interface region of the nearly edge-on 

p ,! M17SW nebula. 

I ■ Methods. We have used the dual color multiple pixel receiver CHAMP+ on APEX telescope to obtain a 5'.3 x 4'. 7 map of the 

O' 7 = 6-^5 and 7 = 7^6 transitions of '"CO , the ''CO 7 = 6 -> 5 line, and the ^Pi -> 370 iim fine-structure transition of [C 

^ I I] in M17 SW. LTE and non-LTE radiative transfer models are used to constrain the ambient conditions. 

C/3 , Results. The warm gas extends up to a distance of ~ 2.2 pc from the M17 SW ridge. The "CO 7 = 6 — » 5 and [C I] 370 yum lines 
have a narrower spatial extent of about 1.3 pc along a strip line at P.A=63°. The structure and distribution of the [C I] 'P? — > '^i 370 

' yum map indicate that its emission arises from the interclump medium with densities of the order of 10' cm"', 

-si I Conclusions. The warmest gas is located along the ridge of the cloud, close to the ionization front. An LTE approximation indicates 

►-^ , that the excitation temperature of the embedded clumps goes up to ~ 120 K. The non-LTE model suggests that the kinetic temperature 

I ■ at four selected positions cannot exceed 230 K in clumps of density n(H2) ~ 5 x 10^ cm"', and that the warm (r^ > 100 K) and dense 

' ("(Ha) > lO** cm"') gas traced by the mid-7 '^CO lines represent just about 2% of the bulk of the molecular gas. The clump volume 

Q,^ I filling factor ranges between 0.04 and 0.1 1 at these positions. 

I Key words, galactic; ISM — galactic: individual: M17 SW — radio lines: galactic — radio fines: ISM — molecules: ''CO , "CO — 

, atoms: [C I] 

^ ■ 

On : 

1. Introduction may also be an important source ofheating in high velocity wing 

sources like, Orion, W5 1 and W49 (Jaffe et al. 1987). 

. !Z| The heating and cooling balance in photon-dominated regions 

, (PDRs) remains an active study ofresearch. The comprehensive Because of its nearly edge-on geometry, and the large 

<J . understanding of PDRs requires observations of large areas close ^"1°^"^ °^ observational data available in the literature, M17 SW 
a ■ to radiation sources, and of a wide wavelength range covering ""^ °^ *e best Galactic regions to study the entire structure of 
■ " " various emissions of atoms, molecules, and grains. In particular, P^Rs from the exciting sources to the ionization front, and the 
mid-7 CO lines have been detected in almost all known massive succession (or not) of Hz, [C I] and CO emissions, as predicted 
Galactic star forming regions (e.g. Orion Nebula, W51, Cepheus ^y PDR models (Icke, Gadey, & Israel 1980; Felli, Churchwell, 
A, NGC 2024). This indicates that warm (Tk > 50 K) and dense & ^assi 1984; Meixner et al. 1992; Meijerink & Spaans 2005). 
(n(H2) > 10' cm"3) gas is common, and probably of importance M17 SW is also one of the few star forming regions for which 
in most OB star forming regions. The mid-7 CO lines detected *e magnetic field strength can be measured m the PDR inter- 
in regions like, e.g. M17, Cepheus A and W51, have relatively f^'^e, and where the structure of the neutral and molecular gas 
narrow line widths of 5-10 km s"' , although not as narrow as the ^eems to be dominated by magnetic pressure rather than by gas 
Une widths observed in cold quiescent cloud cores. pressure (Pellegrini et al. 2007). 

Observations of the 7 = 6 ^ 5 and 7 = 7 ^ 6 transitions M17 SW is a giant molecular cloud at a distance of 2.2 kpc, 

of '^CO in several massive star forming regions indicate that the illuminated by a highly obscured (A,. > 10 mag) cluster of sev- 

warm emitting gas is confined to narrow (< 1 pc) zones close eral OB stars (among k 100 stars) at about 1 pc to the East (Beetz 

to the ionization front. These observations favor photoelectric et al. 1976; Hanson et al. 1997). It also harbours a number of 

heating of the warm gas by UV radiation fields outside the HII candidate young stellar objects that have recently been found 

regions (e.g. Harris et al. 1987; Graf et al. 1993; Yamamoto et (Povich et al. 2009). Several studies of molecular emission, ex- 

al. 2001; Kramer et al. 2004 and 2008). Nevertheless, shocks citation and fine profiles (e.g. Snell et al. 1984; Martin, Sanders 

& Hills 1984; Stutzki & Gusten 1990) from the M17 SW core in- 
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homogeneous. Emission of [C I] and [C II] was detected more 
than a parsec into the molecular cloud along cuts through the 
interface region (Keene et al. 1985; Genzel et al. 1988; Stutzki 
et al. 1988). These results, as well as those found in other star 
forming regions like S106, the Orion Molecular Cloud, and the 
NGC 7023 Nebula (e.g. Gerin & Phillips 1998; Yamamoto et al. 
2001; Schneider et al. 2002, 2003; Mookerjea et al. 2003) do 
not agree with the atomic and molecular stratification predicted 
by standard steady-state PDR models. However, the extended [C 
I] ^Pi and '^^CO 7 = 2^1 emission in S140 have been 

successfully explained by a stationary, but clumpy, PDR model 
(Spaans 1996; Spaans & van Dishoeck 1997). Hence, the lack of 
stratification in [C I] , [C II] and CO is a result that can be ex- 
pected for inhomogeneous clouds, where each clump acts as an 
individual PDR. On the other hand, a partial face-on illumination 
of the molecular clouds would also suppress stratification. 

Based on analysis of low-J lines of '^CO , '^CO and 
CH3CCH data, the temperature towards the M17 SW cloud core 
has been estimated as 50-60 K, whereas the mean cloud tem- 
perature has been found to be about 30-35 K (e.g. Giisten & 
Fiebig 1988; Bergin et al. 1994; Wilson et al. 1999; Howe et al. 
2000; Snell et al. 2000). Temperatures of ~ 275 K has been esti- 
mated from NH3 observations (Giisten & Fiebig 1988) towards 
the VLA continuum arc, which agree with estimates from highly 
excited '^CO transitions (HaiTis et al. 1987). Multitransition CS 
and HC3N observations indicates that the density at the core re- 
gion of M17 SW is about 6 x 10^ cm^^ (e.g., Snell et al. 1984; 
Wang et al. 1993; Bergin, Snell & Goldsmith 1996). While den- 
sities up to 3 X 10^ cm"^ have been estimated towards the north 
rim with multitransition observations of NH3, which indicates 
that ammonia is coexistent with high density material traced in 
CS and HCN (Gusten & Fiebig 1988). The UV radiation field 
Go has been estimated to be of the order of 10"* in units of the 
ambient interstellar radiation field ( 1 .2 x 1 0"^ ergs s" ' cm" ' sr" ' , 
Habing 1968; Meixneret al. 1992). 

However, most of the millimeter-wave molecular observa- 
tions in M17 SW are sensitive only to low temperatures (< 100 
K), and the few available data of mid-y CO and [C I] lines (con- 
sisting mostly of cuts across the ionization front and observa- 
tions at few selected positions) are limited in spatial resolution 
and extent (e.g. Hams et al. 1987; Stutzki et al. 1988; Genzel et 
al. 1988; Stutzki & Gusten 1990; Meixner et al. 1992; Graf et 
al. 1993; Howe et al. 2000). Therefore, in this work we present 
maps (~ 3 X 3 pc^) of mid-7 molecular ('^CO and '^CO ) and 
atomic ([C I] ) gas, with excellent high resolution (9.4" - 7.7"), 
which advances existing work in Ml 7 SW. 

The observations were done with CHAMP^ (Carbon 
Heterodyne Array of the MPIfR) on the Atacama Pathfinder 
Experiment (APEXQ) (Gusten et al. 2006). The multiple pix- 
els at two submm frequencies of CHAMP+, allow for efficient 
mapping of ~arcmin regions, and provide the ability to observe 
simultaneously the emission from the 7 = 6 — > 5 and / = 7 — > 6 
rotational transitions of '^CO at 691.473 GHz and 806.652 GHz, 
respectively. We also observed the 7 = 6 — » 5 transition of 
'^CO at 661.067 GHz and the ^Pi 370 fim (hereafter: 

2^1) fine-structure transition of [C I] at 809.342 GHz. 

Since the gas phase cools mainly via the atomic fine struc- 
ture lines of [O I] , [C II] , [C I] , and the rotational CO lines 
(e.g. Kaufman et al. 1999, Meijerink & Spaans 2005), these car- 



' This publication is based on data acquired with thie Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



bon bearing species presented here are very important coolants 
in the interstellar medium (ISM) of a variety of sources in the 
Universe, from Galactic star forming regions, the Milky Way as 
a galaxy, and external galaxies up to high redshifts (e.g. Fixsen 
et al. 1999; Weiss et al. 2003; Ki-amer et al. 2005; Bayet et al. 
2006; Jakob etal. 2007). 

The case of M17 SW can be considered as a proxy for extra 
galactic star forming regions. M17 SW is not special, nor does 
it need to be, compared to other massive star-forming regions 
like Orion, W49, Cepheus A, or W51. Still, it does allow feed- 
back effects, expected to be important for starburst and active 
galaxies, to be studied in great spatial detail. Comparison of the 
local line ratios to extra-galactic regions can then shed light on 
the properties of massive star forming regions that drive the en- 
ergetics of active galaxies. Our results will be of great use for 
future high resolution observations, since molecular clouds of 
the size of the maps we present will be resolved by ALMA at 
the distance (~ 14 Mpc) of galaxies like NGC 1068. 

The main purpose of this work is to explore the actual spatial 
distribution of the mid-y '^CO and [C I] lines in M17 SW, and 
to test the ambient conditions of the warm gas. A simple LTE 
model based on the ratio between the '^CO and '^CO 7 = 6 — > 5 
lines is used to probe the temperature of the warm (Tk ~ 100 
K) and dense («h > 10^ cm"^) molecular gas. Then a non-LTE 
model is used to test the ambient conditions at four selected po- 
sitions. In a follow up work we will present an elaborate model 
of these high resolution data. 

The most frequent references to Stutzki et al. (1988), Stutzki 
& Gusten (1990) and Meixner et al. (1992) will be referred to as 
S88, SG90 and M92, respectively. The organization of this arti- 
cle is as follows. In §2 we describe the observations. The maps 
of the four lines observed are presented in §3. The modelling and 
analysis of the ambient conditions are presented in §4. And the 
conclusions and final remarks are presented in §5. 



2. Observations 

We have used the dual color heterodyne array receiver CHAMP+ 
(Kasemann et al. 2006; Gusten et al. 2008), providing 2x7 pix- 
els, on the APEX telescope during July 2008, to map simulta- 
neously the 7 = 6 — » 5 and J - 7 6 lines of '^CO , 
and - in a second coverage - the i^CO 7 = 6 ^ 5 and [C 
I] 7 = 2 ^ 1. We observed a region of about 5'. 3 x 4'.7 (3.4 
pc X 3.0 pc) in on-the-fly (OTF) slews in R.A. (~ 320 arc- 
sec long), subsequent scans spaced by 4" in Declination. The 
observations were done in total power mode (nodding the an- 
tenna prior to each OTF slew to a reference position 180" east 
of the SAO star 161357. The latter is used as reference through- 
out the paper, with R.A(J2000)=18:20:27.64 and Dec(J2000)=- 
16:12:00.90. We used Sgr B2(N) as reference for continuum 
poiting. Calibration measurements were performed regularly, 
every ~10 min, with a cold liquid nitrogen (LN2) load and an 
ambient temperature load. The data were processed with the 
APEX real-time calibration software (Muders et al. 2006), as- 
suming an image sideband suppression of 10 dB. 

We used the Fast Fourier Transform Spectrometer (EFTS) as 
backend with a fixed bandwidth of 1.5 GHz and 1024 channels. 
We used the two IF groups of the EFTS with an offset of +460 
MHz between them. The spectral resolution was smoothed to 
about 1 kms"', while the line widths are between 4 kms"' and 
9 kms so they are well resolved. The on source integration 
time per dump and pixel was 1 second only. However, oversam- 
pling with 4" spacing, all the seven pixels of CHAMP+ covered 
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Fig. 2. Color map of the 21 cm continuum emission (Jy beam ') in 
M17 SW with 10" X 7" resolution by Brogan & Troland (2001). The 
white contour lines correspond to the '^CO 7 = 6 — » 5 with a peak 
emissions of 830 K kms"'. While the green contour lines correspond 
to the [CI] J = 2 —> 1 with a peak emissions of 260 K kms"'. The 
countour levels (from thin to thick) are the 25%, 50%, 75% and 90% 
of the peak emission. The red stars indicates the O and B ionizing stars 
(Beetz et al. 1976; Hanson et al. 1997). The reference position (Aa = 0, 
A(5 = 0) is the same as in FigE These '-CO and [C I] maps have a 
slightly lower integrated temperature than in Fig.[T]because they were 
convolved with a 20" beam to smooth the contour lines. 
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Fig. 1. Top - Color map of the integrated temperature of '^CO } = 

6 ^ 5 in M17 SW. The contour lines correspond to the '^CO J = 

7 — » 6, which has a peak emission of 925 K kms"'. The countour 
levels are the 25%, 50%, 75% and 90% of the peak emission. Dashed 
lines correspond to the E-W and NE-SW strip lines at P.A=90° and 
P.A=63°, respectively. Bottom - Color map of the integrated temper- 
ature of '^CO y = 6 ^ 5 and the contour levels (as described be- 
fore) of [C I] / = 2 — > 1 with a peak emission of 282 K kms"'. 
The filled triangle and square mark selected positions where ambi- 
ent conditions are estimated from. The reference position (Aa = 0, 
AS = 0), marked with a cross, corresponds to the SAO star 161357 
at R.A(J2000)= 18:20:27.6483 and Dec(J2000)=-16: 12:00.9077. 



a given grid position at least once. So, after adding all the sub- 
scans from both IF channels, and after convolving the maps with 
the corresponding beam size, the total integration time in the 
central 5'x4' region of the maps varies between about 50 and 80 
seconds per grid cell. 

The SSB system temperatures are typically about 2000 K 
and 6000 K for the low and high frequency bands, respectively. 
The spatial resolution varies between 9.4" for '^CO 7 = 6^5 
transition in the low frequency band (at 661 GHz - the nomi- 
nal beam at 691 GHz is 8.4") and 7.7" for the high frequency 
band (809 GHz). All data in the paper were converted to line 
brightness temperature, Tb = rjf x T^/rjc, using a forward effi- 
ciency (rif) of 0.95 and beam coupling efficiencies (rjc) of 0.45 
and 0.43 (at 661 GHz and 809 GHz, respectively) as determined 



towards JupiteiQ (Giisten et al. 2008). We assumed brightness 
temperatures of 150 K (at 660 GHz) and 145 (at 815 GHz) for 
Jupiter (Griffin et al. 1986). This coupling efficiency was cho- 
sen because in velocity-space (velocity channels) the size of the 
M17 clumps is Jupiter-like, which had a size ~ 38.7" by the 
time of the observations. The calibrated data were reduced with 
the GILDAS0 package CLASS90. 



3. Results 

3.1. Integrated line temperature maps 

Figure [T] shows the maps of the temperature, integrated between 
5 kms"' and 35 kms"', of '^CO 7 = 6^5 (top) with the 
contour lines of '^CO 7 = 7-^6, and the velocity integrated 
temperature of '-^CO 7 = 6 — > 5 (bottom) with the contour lines 
coiTesponding to [C I] 7 = 2 ^ 1 . All the maps were convolved 
to the largest beam size (9.4") of the '^CO 7 = 6^5 line, 
obtaining a grid size of about 4.7" x 4.7". The peak integrated 
temperature of the '^CO 7 = 6^5 and 7 = 7 — > 6 lines are 
852 K km s"' and 925 K km s"' respectively. These lines follow 
a similar spatial distribution. The peak integrated temperature 
of i^CO 7 = 6 ^ 5 and [C I] 7 = 2 -» 1 are 420 K and 
282 K. respectively, and the peak of [C I] is shifted towards the 
inner side of the interface region at about 0.55 pc (~ 50"). The 
ionization front traced by the high resolution (10" x 7") map of 
the 21 cm continuum emission (Brogan & Troland 2001), as well 
as the ionizing stars identified by Beetz et al. (1976) and Hanson 
et al. (1997), are shown in Figure |2l with '^CO 7 = 6^5 
(white contour lines) and [C I] 7 = 2 — » 1 (green contour lines) 
overlaid. The transition between the hot (T^ > 300 K) atomic gas 



' http://www.mpifrde/div/submmtech/heterodyne/champplus/ 
champmain.html 
' http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 3. Top panel - Strip lines of the velocity integrated intensities of 
'^CO 7=7-6 (black), '^CO 7=6-5 (red), '^CO 7=2-1 (green) (adapted 
from S88), "CO 7=6-5 (blue) and [C I] 7=2-1 (cyan) at ^6 = 30" 
(P. A. 90°) across the ionization front of M17 SW. Bottom panel - Strip 
lines at PA. 63° (A(5 = Aa/1). The X-axis corresponds to the actual 
offset in R.A. of the maps shown in Figure [T] Hence, Aa = 0" is the 
R.A. of the reference illuminating star SAO 161357. 



and the warm (Tk > 100 K) molecular gas can be seen thanks to 
the almost edge-on geometry of Ml 7 SW. 

The top panel of Figure [3] shows the variation of the inte- 
grated temperature of all the lines, across the ionization front 
(strip line at P.A=90° in Fig. [TJ- Due to the limited S/N the 
'^CO 7 = 7^6 and [C I] 7 = 2 — > 1 strip lines have been 
smoothed spatially with respect to the strip direction. The [C 
I] 7 = 2 — » 1 line starts peaking up at about 0.1 pc (~ 10") 
after the molecular lines and presents a smooth transition to- 
wards the inner part of the cloud, forming a plateau up to about 
Aa = -100", from where it increases its emission until the peak 
reached at about Aa = -120". The peak of [C I] correlates with 
a secondary peak seen in '^CO . However, the main peak emis- 
sion of the latter correlates with the peak of the '^CO lines along 
this strip line. 

The strip line at P.A.=63° {bottom panel of Figure[3]) can be 
compared with Fig. 5 in M92, and Fig. 2 in S88. At this position 
angle, there is not a marked plateau in the [C I] emission, and 



the peak of the [C I] line is closer to the peaks of the '^CO and 
'^CO lines. The dip in '^CO 7 = 2 ^ 1 at about Aa = -150" is 
an artifact. The integrated temperature of the '^CO 7 = 2^1, 
7 = 6 — > 5 and J - 7 6 lines have comparable strength deep 
(Aa > -160") into the M17 SW complex. 



3.2. The complex internal structure of M17 SW 

Figure |4] shows the spectra at selected positions along the NE- 
SW strip line at P. A. 63°. The main beam temperature of the 
spectra is shifted 70 K at each offset position. This set of spectra 
can be compared with the '^CO and C'**0 7 = 2 — > 1 spectra 
along the same strip line of Fig. 8 in S88. The warm gas {Tfc > 50 
K), traced by the mid-7 '^CO lines, is as extended as the cold 
gas (Tk < 50 K) traced by the '^CO 7 = 2^1 line deeper 
into the cloud. On the other hand, the '^CO 7 = 6^5 and [C 
I] 7 = 2 — > 1 lines are strongly detected in a narrower spatial 
extent of about 1.3 pc, similar to the extent of the C'^O 7 = 2 — > 
1 emission. 

Multilevel molecular line observations in CS, '^CO , 
'^CO and C'^O, and in several fine structure lines ([C I] , [C II] , 
[Si II], [O I]) indicate that M17 SW consists of numerous high 
density clumps («(H2) > 10"* cm"^) where the [O I], [Si II] and 
mid-7 CO lines emanate from. This dense gas is found within 
relatively warm (~ 50 K) and less dense (n(H2) ~ 3 x lO^* cm"-') 
molecular gas (interclump medium) which in turn is surrounded 
by a diffuse halo (n(H2) ~ 300 cm"^) where the very extended 
[C I] and [C II] emission emerge from (Snell et al. 1984, 1986; 
Evans et al. 1987; S88; SG90; M92). 

From J8o observations in M17 SW a beam-averaged (13") 
column density of ~ 8 x 10^^ cm"^ has been estimated for the 
cloud core, and masses in the range ~ 10 - 2000 Mq for the CO 
clumps (SG90). A comparable mass range (~ 10-120 Mq) was 
lately estimated from submillimeter continuum observations in 
the northern part of M17 (Reid & Wilson 2006). Although the 
region mapped by Reid & Wilson (2006) adjoins, but does not 
overlap with M17 SW. 

Figure |5] shows representative velocity channel maps of the 
'^CO 7 = 6^5 {top left) and J - 7 ^ 6 {top right) lines in 
M17 SW. These are the main-beam brightness temperature av- 
eraged over 2 and 3 velocity channels between 18.2 kms"' and 
19.9 km s '. These are similar velocity channels shown in Fig. 3 
by SG90. The fact that the C"*0 7 = 2 ^ 1 line trace colder 
{Tk < 50 K) and less dense («// ~ 3 x 10^ cm"-') gas than the 
'^CO lines is reflected in the different velocity integrated and 
channel maps of these lines. In theory the critical densities (at 
Tk = 100 K) of the '^CO 7 = 6^5 and 7 = 7^6 lines 
are n„;, ~ 2.7 x 10^ cm"-' and ncra ~ 4.4 x 10^ cm"-', respec- 
tively, which corresponds to a factor ~ 1 .6 difference. However, 
this difference does not translates directly into different clumpy- 
ness. This is reflected in the similar clumpy structure seen in the 
channel maps of these mid-7 '^CO lines. 

Even though the critical density of the '^CO 7 = 6 — » 5 
Une is similar to that of the '^CO («„•/, ~ 2.4 x 10^ cm"^) the 
South-East region of its channel map {bottom left) differs from 
that seen with the '^CO lines. This could be due to a change in 
the temperature of the gas, or to a variation in the '^CO column 
density in that region. Since '^CO is much more optically thin 
than '^CO (abundance ratio of about 50 - 70) this difference in 
the map can be expected. In Sections 4.2 and 4.3 we discuss 
about the optical depths. 

On the other hand, the [C I] 7 = 2 — » 1 channel map {bot- 
tom right) shows a completely different structure and distribu- 
tion than the '^CO and the isotope lines. Since the critical den- 
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Fig. 4. Selected spectra of '^CO J = 6 5, '^CO / = 7 ^ 6, "CO 7 = 6 ^ 5 and [C I] 7 = 2^1 along the NE-SW strip line (see Fig.[B. 
The position in arcsecs are the offsets with respect to the reference position at R.A(J2000)= 18:20:27.6483 and Dec(J2000)=-16:12:00.9077. The 
spectra are the average spectra within +2" of the indicated offset positions, and smoothed to 1 kms"' velocity resolution. 



sity of this line is about 2.8 x 10'' cm ^, its emission is likely 
emerging partly from the interclump medium mentioned above. 

4. Discussion 

4.1. Self-absorption in the mid- J ^^CO lines? 

The complex structure of the '^CO J - 1 -^0, 7 = 2— > 1, 
and 7 = 3 — » 2 line profiles has been attributed to strong self- 
absorption effects (e.g. Rainey et al. 1987; Stutzki et al.l988). 
Martin, Sanders & Hills (1984) also reported a flat topped spec- 
tra of '^CO 7 = 3^2, attributed to self-absorption or saturation 
at velocities near the line center, and give details about the effects 
of macroturbulent clumpy medium in line profiles. 

A double peaked structure in the '^CO 7 = 1 — > line was 
also reported by Lada (1976). Rainey et al. considered that this 
double peaked structure in '^CO suggest that either this line is 
optically thick, or that the double peaked structure is due to more 
than one cloud component. The latter is the interpretation fa- 
vored by Rainey et al. in view of the available data at that time. 

Phillips et al. (1981) presented a self-absorption LTE model 
that considers a '-CO cloud of uniform temperature Tk in front 



of a hot background source of temperature Tbg, at the same cen- 
tral velocity. The velocity dispersion of the background cloud is 
considered to be larger compared to that of the foreground cloud, 
so the self-absorption effect is seen mostly at the line center. 
This model indicates that, depending on the total column den- 
sity of '^CO , the self-absorption effect will be stronger in the 
7 = 2 — > 1 and 7 = 3 — > 2 lines than in the 7=1^0 line, with 
decreasing intensity as the transition number 7 increases. This 
is indeed observed in Fig. 12 of S88, for the '^CO 7 = 4 -> 3, 
7 = 3^2, 7 = 2— > 1, and 7 = 1 — > lines. 

We reproduced the model by Phillips et al. including the 
higher-7 lines of '^CO . The top panel of Figure |6] shows the 
model with the same background and foreground temperatures 
used by Phillips et al.. This model implies that, for a background 
temperature Ti,g = 64 K and a foreground kinetic temperature 
Tk - 15 K, the lower-7 lines (J=l,2,3,4) of the background 
cloud start showing self absorption at the line center for lower 
column densities (N/AV = lO'"' - 10'^ cm"^ km s"'). Instead, 
the higher-7 lines (J=5,6,7) need larger columns (N/AV = 10'^ - 
10'^ cm"^ km s"') in order to be affected by self-absorption. For 
a velocity dispersion AV - 5 kms the upper limits of these 
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Fig. 5. Top panel - Channel maps of the main-beam brightness tem- 
perature of '^CO y = 6 — > 5 {left) and 7 = 7 — > 6 (right) averaged 
over 2 and 3 velocity channels, respectively, between 18.2-19.9 km s"'. 
The contour lines are as described in Fig.l, with peak integrated line 
temperatures of 54.5 Kkm"' s and 55.7 Kkm"' s for the 7 = 6 — > 5 
and 7 = 7 — > 6 lines, respectively. Bottom panel - Channel maps of 
"CO y = 6 -> 5 (left) and [C I] 7 = 2 ^ 1 {right), averaged over 
2 velocity channels between 18.4-19.7 kms"'. Contours are as in top, 
and the peak integrated line temperatures are 37.0 Kkm"' s and 23.9 
Kkm ' s . 



^^CO columns would correspond to extinctions A,, of ~ 0. 1 mag 
and ~ 10 mag, respectively. 

The bottom panel of of Figure |6] shows the model for a 
background temperature T^g - 150 K and a foreground tem- 
perature = 30 K (from S88). In this case the lower-/ lines 
show self-absorption at the same range of columns as before, 
while the higher-/ lines start showing self-absorption at a nar- 
rower range of columns (A^/AV^IO'^ - lO'*- cm^^ km s"'). A 
remarkable characteristic of these models (top and bottom pan- 
els of Figure |6]l is that all the J lines are expected to be strongly 
self-absorbed at columns larger than 10^^ cm"^ km s which is 
similar to the column density estimated by S88. Another char- 
acteristic is that the '^CO emission of the higher-7 lines are also 
expected to decrease with the transition number J, and be weaker 
than the low-7 lines. However, the '^CO J - 7 —> 6 line seems 
to break this rule, as can be seen in fig. 12 of S88. The high peak 
temperature observed in the '^CO 7 = 7 — > 6 line is missing 
in the lower-/ lines. Even considering a calibration uncertainty 
of 20%, the '^CO / = 7 — > 6 line (observed at offset position 
(-100", 0"), bottom panel of Fig. 12 in S88) will be as strong 
as the 7 = 4 — > 3 line (at least at the peak intensity) but still 
stronger than the 7 = 2 — > 1 line. 

On the other hand, the '^CO 7 = 7 — > 6 line seems to be 
asymmetric, with a left shoulder weaker than the right shoulder, 
which may be due to self-absorption produced by a colder fore- 
ground cloud with slightly lower center velocity than the warmer 
clump traced by the J - 7 ^ 6 line. However, that weaker left 
shoulder of the mid-7 line is still brighter than the correspond- 
ing shoulder of the lower-7 lines, in most of the velocity range 
and in both positions (-100", 0") and (-60", -30") - assuming 
low (< 10%) uncertainty in the calibration of the data. This is not 
what would be expected in the self-absorption scenario proposed 
by Phillips etal. (1981). 




log^^(N/AV) [cm ^km s '] 




Fig. 6. Top panel - Expected brightness temperature at the center of the 
'^CO lines for a warm background cloud with temperature Ti,g = 64 
K and a colder foreground absorbing cloud with temperature T); = 15 
K. Bottom panel - Same as in top, but for a background temperature 
Ttg = 100 K and a foreground temperature of Tj. = 30 K. 



Figure m shows that the '^CO 7 = 6-^5 line has a sim- 
ilar asymmetry as the '^CO lines, at positions (-45", -23") 
and (-60", -30"). But it shows only one component at the 
other positions. This difference may be related to a gradient 
in the temperature (or total column density) of the foreground 
cloud that produces self-absorption in the first two positions, 
but not in the others. Instead, [CI] 7 = 2 — > 1 shows similar 
asymmetry as '^CO at positions (-45", -23"), (-150", -75") 
and (-180", -90"), and an opposite asymmetry at position 
(-75", -38"). Given that there is no strong evidence for self- 
absorption in the '^CO lines, nor in the [C I] lines, and that 
the '^CO lines are mostly optically thin, it is unlikely that the 
observed asymmetries of the '^CO and [C I] lines are produced 
by self-absorption. Hence, we agree with Rainey et al. (1987) in 
that this complex structure is more likely due to more than one 
kinematical component along the line of sight. And this could 
also be the case for the mid-7 '^CO lines. 

Therefore, the observational facts and the models suggest 
that the self absorption effect, if present, should have little im- 
pact on the mid-7 lines, and a few cloud components at dif- 
ferent central velocities could also explain the complex struc- 
ture of the line profiles. The asymmetry of the profiles suggests 
that self-absorption affects mostly one wing of the line profile, 
while the peak temperatures seems to be the least affected veloc- 
ity channel in the mid-7 lines. Hence, in the following sections 
we test the ambient conditions of the warm gas based on the ra- 
tios between the peak main-beam temperatures of the '^CO and 
i^CO 7 = 6^5 and 7 = 7->6 lines. 
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4.2. Optical depth and excitation temperature (LIE) 

Since we have the maps of '^CO and '^'CO 7 = 6 — > 5 lines, we 
can estimate the optical depth and the excitation temperature of 
these lines, assuming local thermal equilibrium (LTE), from the 
ratio between their peak main beam temperature T„,h observed 
between 5 km s ' and 35 km s ' velocity channels. This will pro- 
vide at least a lower limit for the kinetic temperature in M17 SW. 
Then we will estimate the ambient conditions at two selected po- 
sitions based on a non-LTE model of the ratio between the peak 
T„jh temperatures of the '^CO 7 = 6 — > 5 and J = 7 ^ 6 lines 
(hereafter referred as '^CO ^ line ratio). The temperature and 
densities obtained in this way will be compared to those values 
estimated in previous work. 

In LTE the radiation temperature can be approximated (e.g. 
Kutner 1984; Bergin et al. 1994) by the expression: 

r« = [7v(7'«)-7bg][l-e-^1, (1) 

where the term JyiT) is the Planck's function evaluated at fre- 
quency V and temperature T, and multiplied by the factor to 
obtain the intensity in K. So it is defined as: 



hv, 



'A 



■"/AT _ I ■ 



(2) 



We use the full Jv(T) function since the Rayleigh- Jeans (R- 
J) approximation (commonly applied when hv «: kT) does not 
hold for the high frequency lines studied in this work. For the 
R-J approximation to be true, we require T » 300 K, which is 
much higher temperature than what we expect to trace with our 
observations. 

The background radiation 7bg is a composite between the 
cosmic microwave background radiation (CMB), as a blackbody 
function at 2.73 K, and the diluted infrared radiation remitted by 
dust. That is: 



7bg = 7v(2.73) + TdJyiTd), 



(3) 



where is the effective optical depth of the warm surface layer, 
adopted from Hollenbach et al. (1991), and it is defined as 
Td - Ti(%m(100yum//l). For M17 SW we adopted an emission 
optical depth at 100 yum of rioo/zm = 0.106 and the average dust 
temperature Td - 50 K from M92. We tried both, with and with- 
out the dust contribution to the background radiation, and we 
found that the contribution of the radiation by dust continuum 
emission is negligible at frequencies of the order of 690 GHz 
and 810 GHz. Nevertheless, all the following analysis include 
the dust contribution for completeness. 

For extended (resolved) sources like the clumps in M17 SW, 
the radiation temperature is well estimated by the observed main 
beam brightness temperature T„,i,. Hence, we use that quan- 
tity in the following analysis. From the LTE approximation we 
can assume that the excitation temperatures T^^ of '^CO and 
'^CO 7 = 6 — > 5 are the same, although the terms JyiJex) are 
not exactly the same because of the slightly (~ 4%) different fre- 
quencies of the '^CO and '^CO lines. So, from equation (1) the 
ratio between '^CO and '^CO can be approximated as: 



T„,hC^COJ ^ 6-5) 



I _ g-T('-CO 7=6-5) 

6-5) ~ 1 _ e-T("co 7=6-5) ' 



(4) 



Following the work by Wilson et al. (1999), we adopted 
a constant ['^CO /'^CO ] abundance ratio of 50 for M17 
SW, which is approximately the value measured at a similar 



100 



-100 




Fig. 7. LTE approximation of the optical depth {Top) of the '^CO J = 
6 — > 5 and the excitation temperature (Bottom) of the "CO and 
'^CO 7 = 6^5 lines in M17 SW. Contours ai'e the 10%, 25%, 50%, 
75% and 90% of the peak value, which is 1.9 for tO^CO 6-5) and 120 
KforJ,,. 



Galactic radius towards the W51 region (Langer & Penzias 
1990). Assuming that the optical depth is proportional to the to- 
tal column density of the molecules and, hence, to the abundance 
ratio between them, we can estimate that t('-CO) ~ 50r('-'CO). 
The '^CO line is usually optically thin, so r('^CO) could be 
taken out of the exponential in equation (3), and estimated di- 
rectly. However, we do not really know if this holds for the en- 
tire Ml 7 SW region, so we do not apply further approximations 
and we solved equation (3) for t('^CO) with a numerical method 
(Ne wton-Raphson) . 

The top panel in Figure [T] shows the t('^CO) map. The 
'^CO line is optically thin in most of the region, with some 
optically thick spots (e.g., Aa = -30, A5 = -110). Knowing 
t('''CO) we can estimate r„ from equation (2) using either 
tracer, considering that the Tg^ estimated using '^CO is just 
~ 0.6% higher than that estimated using '^CO . The T^x map is 
shown in the bottom panel of Figure [T] This map indicates that 
the warmest gas is located along the ridge of the cloud, close 
to the ionization front. The temperature in this region ranges 
between 40 and 120 K, and the peak temperature is located at 
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around (Aa - -60, A6 - 10). If we consider only the gas with 
temperatures > 80 K, the warm gas would be confined to a zone 
of about 40" (~ 0.44 pc) next to, and along, the Hll region, 
which agrees with previous results found by Graf et al. (1993). 
If the gas were thermalized, then this can be the actual map of 
the kinetic temperature of the gas. Otherwise, the map can 
be considered as a lower limit of Tk- Since in velocity space 
the clumps cover the whole beam, this would imply that the 
^^CO and '^CO molecules are subthermal in the 7 = 6 — > 5 
transition. That is, the density of the gas and the column density 
of ^^CO and '^CO may be insufficient to thermalize these transi- 
tions. A more detailed analysis is presented in the next section. 

4.3. Ambient condition at selected positions (non-LTE) 

Figure [8] shows the spectra of all the observed lines extracted at 
two diff'erent positions in the map. The top panel shows the spec- 
tra observed at position A (Ac = -70", A(5 — +2)2"), close to the 
peak emission of the '^CO lines. All the lines show a double 
component structure with the secondary component peaking at 
~ 25 kms"'. The middle top panel shows the spectra at position 
B (Aa - -70", A5 = -82"), where the velocity integrated tem- 
perature corresponds to about 50% of the peak emission. Here 
only the '^CO 7 = 7 — > 6 line seems to have a deep at the 
line center However, because the low S/N in the high frequency 
band, this dip may be likely due to noise. The middle bottom 
panel shows the spectra at position C (Aa - -60", A6 - -30"), 
which corresponds to the peak of the NE-S W strip scan reported 
in S88 and Graf et al. (1993), with beams of 40" and 8", re- 
spectively. The bottom panel shows the spectra at position D 
(Aa = -100", A5 - 0"), which is close to the continuum far- 
IR peak, also reported in S88. Since we do not have dedicated 
observations at these positions, we extracted the spectra from the 
nearest pixels in our maps, convolved to the largest beam (9.4") 
of the '^CO y = 6 — » 5 line. So the spectra shown in Figure [8] 
are the convolved spectra centered within +1" of the indicated 
coordinates. This is justified because we have oversampled data. 

Table[T]shows the Gaussian fits of the spectra obtained at the 
four selected positions. Two Gaussian components were needed 
to fit the lines, except at position B, where only one component 
was used. The main components of the '^CO lines have a line 
width that is about 8-9 km s"' at position A, while the '^CO has 
a line width of about 3 km s ' narrower. The [C I] line is the nar- 
rowest line, with a line width of ~ 4 kms '. At position B, the 
^^CO lines are the widest of the four lines with about 8 kms ' 
and the '^CO and [C I] lines have about half the line width of 
the '^CO lines. At position C and D the Gaussian parameters of 
the ^^CO 7 = 7 — > 6 presented uncertainties of ~ 50% when 
let free in the fitting. However, because the line shape of the 
'^CO J = 1 6 and 7 = 6 — > 5 transitions are very similar, 
we set the line width of the 7 = 7 — > 6 transition to the value 
found for the 7 = 6 — » 5 line. The line width of the main com- 
ponents of the '^CO lines at position C and D are ~ 6kms 
That is, about 2 kms ' narrower than the lines observed at po- 
sitions A and B. This difiference can be due to a higher optical 
depth towards the latter positions, or to the contribution of few 
fast-moving cloudlets (Martin, Sanders & Hills 1984; Graf et al. 
1993). 

The '^CO 1^ line ratio between the peak main beam tem- 
peratures T„b obtained from the Gaussian fit of the main com- 
ponents is 1.02 + 0.05 at position A, 0.95 + 0.05 at position B, 
1.00 + 0.05 at position C, and 0.99 + 0.07 at position D. From 
line ratios we can estimate the ambient conditions for these par- 
ticular positions. We have used the non-LTE radiative transfer 



code RADEXQ (Van der Tak et al. |2007l l for estimating the av- 
erage ambient conditions (kinetic temperature, density and col- 
umn density) of the molecular gas. We assumed collisional ex- 
citation by molecular hydrogen. We also assumed an homoge- 
neous spherical symmetry in the clumps for the escape probabil- 
ity formalism. The collision rates between '^CO and ortho- and 
para-H2 are taken from Wernli et al. (2006), and can be found 
in the LAMDA database (Schoier et al. 2005). As in the LTE 
case, we used the cosmic microwave background radiation at 
2.73 K, and we also tested the non-LTE model with and without 
the infrared radiation remitted by dust (eq.3) as the background 
source. It was also found that dust continuum emission produces 
a negligible effect in the non-LTE model at the frequencies of 
the '^CO 7 = 6 — » 5 and J - 7 ^ 6 lines. We explored molec- 
ular hydrogen densities between 10"* cm"^ and 10^ cm"^, tem- 
peratures between 5 K and 500 K, and '^CO column densities 
between 10'" cm"^ and 10"* cm"^. 

Figure |9] shows the possible ambient conditions required to 
reproduce the '^CO line ratios, and the peak T^h of the 
'^CO 7 = 6^5 line observed at position A (top panel) and B 
{bottom panel). A wide range of temperatures (100 - 450 K) and 
densities (> 3 x 10'* cm"^) are possible solutions for a '^CO col- 
umn density per line width A^('2cO)/Ay ~ 5x10'*' cm km s. 

Figure [10] shows the possible ambient conditions estimated 
for position C (top panel) and D (bottom panel). The combina- 
tions of temperatures and densities required to reproduce the line 
ratios and peak temperatures are similar to those found for posi- 
tion A and B. Although, the range of possible temperatures (for a 
given density) at position D is larger than at the other positions. 
The column densities differ due to the different line strengths 
observed at the four positions (Table[TJ. 

In order to constrain the range of solutions we can adopt the 
average 5x10^ cm"-^ density estimated by M92, which is also 
similar to the mean density of the clumps estimated by SG90. 
This is a sensitive assumption for a collision dominated sce- 
nario since this density is larger than the critical density of both 
'^CO hnes for Tk > 20 K. However, at this density (5x10^ cm"^) 
the temperature cannot be higher than 230 K in order to repro- 
duce the line ratio and the peak T^i, of the '^CO 7 = 6 — > 5 
line observed at position A. And it cannot be higher than 150 K 
at position B. At position C the limit is about 220 K, and at po- 
sition D it is about 200 K. These are lower kinetic temperatures 
than the 1000 K estimated for the dense clumps in the three com- 
ponent model proposed by M92. Our upper limits for the kinetic 
temperature agree with the results reported in previous work (e.g 
Harris et al. 1987; S88; SG90). From the map of the excitation 
temperature j estimated from the LTE model (Figure |7]i the 
lower limits for Tk would be ~ 110 K and ~ 80 K at position 
A & C and B & D, respectively. These are similar (within 30%) 
to the lowest temperatures obtained with the non-LTE models 
(Figures|9]&[l0ll. 

According to the radiative transfer model temperatures up to 
400 K, and higher, are also possible, but they require densities 
< 10^ cm"^ in order to reproduce the observed line ratios and 
peak temperatures. These densities and temperatures are consis- 
tent with the estimates made based on previous observations of 
the I^CO 7 = 7 ^ 6 and 7 = 14 -> 13 lines (Hams et al. 1987; 
SG90). On the other hand, clumps with densities > 10'' cm"^ 
could also reproduce the observed ratios and peak T„,i, in all the 
positions, at temperatures < 200 K. However, these would be at 
the lower limit of the temperature range estimated in Harris et 
al. (1987) and SG90. The densities and temperatures found for 

http://www.sron.rug.nl/~vdtak/radex/radex_manual.pdf 
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CO J=6-5 
CO J=7-6 
'■'cO J=6-5 
[C I] J=2-1 



(-70, +32) 





100 




CO J=6-5 
CO J=7-6 
'■^CO J=6-5 
[C I] J=2-1 



(-70. -82) 



(-60. -30) 




(-100, 0) 




Table 1. M17 SW line parameters derived from Gaussian fits, at 
four selected positions. 



Velocity (km/s) 



Fig. 8. Top - Spectra of the four lines observed in Ml 7 SW at posi- 
tion A (Aa = -70",A<5 = +32"), close to the peak emission in the 
'^CO and "CO maps. Middle top - Spectra observed at position B 
(Aa = -70", A5 = -82"), where the integrated line temperatures are 
about 50% of the peak emission. Middle bottom - Spectra observed at 
position C (Ao- = -60", A(5 = -30"), the peak of the NE-SW strip scan. 
Bottom - Spectra observed at position D (Ace = -100", A5 = 0"), close 
to the continuum far-IR peak. 



Ml 7 SW are similar to those found in W3 Main (Kramer et al. 
2004), but higher (although compatible) than the kinetic temper- 
atures found in Carina, where lower limits between 30 K and 50 
K were estimated (Kramer et al. 2008). 



Molecule-i 


Imb 


V 


mb 


AK 






[Kkms"'] 


[kms"'] 


[K] 


[kms"'] 


Position A (-70", +32") 


'2CO / 


=6-5 


723.1±7.5 
44.6±3.3 


20.4+0.04 
25.1+0.03 


82.9+1.5 
31.5+2.9 


8.2+0.11 
1.3+0.07 


'-CO J 


=7-6 


825.4±24.1 
56.9±12.8 


20.3+0.13 
24.2+0.15 


84.6+3.6 
26.5+7.3 


9.2+0.28 
2.0+0.32 


'^CO J 


=6-5 


301.2±4.2 
51.3±3.6 


20.2+0.04 
23.9+0.04 


54.7+1.1 
21.2+1.7 


5.2+0.08 
2.3+0.09 


[C I] J= 


2-1 


138.5+11.5 
29.9±10.1 


20.2+0.16 
23.6+0.18 


30.0+3.8 
12.7+5.2 


4.3+0.42 
2.2+0.51 


Position B (-70", -82") 


'-CO J 


=6-5 


537.4±6.9 


19.7+0.05 


65.7+1.3 


7.7+0.11 


'-CO J 


=7-6 


528.5+16.1 


19.5+0.12 


62.2+2.8 


7.9+0.26 


'^CO J 


=6-5 


166.7+1.6 


19.3+0.02 


33.5+0.5 


4.7+0.05 


[C I] J= 


2-1 


169.7+4.2 


19.2+0.06 


34.4+1.4 


4.6+0.14 


Position C (-60", -30") 


'2C0 J 


=6-5 


220.4+3.1 
365.8+4.9 


17.9+0.04 
21.9+0.03 


33.7+0.6 
86.2+1.6 


6.1+0.06 
3.9+0.05 


'-co J 


=7-6 


220.2+16.8 
367.1+16.4 


17.9+0.31 
21.9+0.08 


33.6+2.6 
86.5+4.0 


6.1" 
3.9" 


'^CO J 


=6-5 


145.7+11.5 
120.3+11.3 


18.9+0.17 
21.8+0.05 


29.8+2.7 
36.9+3.6 


4.6+0.21 
3.1+0.09 


[C I] J= 


2-1 


41.6+27.6 
124.5+28.6 


17.5+0.55 
20.7+0.48 


12.1+8.6 
24.9+7.3 


3.2+0.77 
4.5+0.72 


Position D (-100", 0") 


'^CO J 


=6-5 


136.1 + 10.9 
374.4+11.6 


16.8+0.09 
21.8+0.09 


34.6+3.2 
57.9+2.5 


3.7+0.17 
6.1+0.18 


'^CO J 


=7-6 


114.1 + 11.5 
369.7+18.9 


16.9+0.19 
21.6+0.15 


29.1+3.2 
57.2+3.4 


3.7" 
6.1" 


'^CO / 


=6-5 


82.7+4.2 
131.3+4.1 


18.4+0.09 
21.6+0.05 


20.0+1.3 
28.7+0.9 


3.9+0.15 
4.2+0.04 


[C I] J= 


2-1 


95.9+5.2 
120.9+1.5 


18.7+0.17 
21.4+0.06 


17.9+1.5 
24.5+1.3 


5.0+0.30 
4.6+0.24 



The uncertainty of this parameter was larger than 50%, when let 
free in the Gaussian fitting. So we set its value to the one found for 
the corresponding Gaussian component of the '-CO 7=6-5 line. 



4.3.1 . Column densities at selected positions 

The column density per line width A^('^CO)/Ay at position A 
varies over a small range of 4 - 6 X 10'^ cm"^km~' s. If we 
consider an average line width of 8.7 kms"' estimated for the 
'^CO y = 6 — > 5 and J - 1 —> 6 lines (Table [TJ and the average 
A^(CO)/Ay = 5 X 10'^ cm"- km"' s, then we have a total column 
density of NC^ CO) ~4x 10'^ cm^^. 



The model indicates that at position B the CO column 
density per line width would be ~ 4 x lO'^ cm"^km"' s. The 
average line width of the lines at position B is ~ 7.7 kms"' 
(Table [1), which gives a total column density A^('^CO) ^ 3 x 
10'^ cm"^, similar to the column found at position A. At posi- 
tions C and D the average column densities per line width are 
~ 8 X 10'^ cm"^km"' s and ~ 5 x 10"' cm"^km"' s, respec- 



10 



J.P. Perez-Beaupuits et al.: Warm gas in Ml 7 SW 





450 




400 




350 


g 


300 






250 








200 




150 




100 




50 



(A) 



(-70", +32") 



16.6 



16.7 



16.8 



CO 1^ = 1.02 ± 0.05 

6-»5 



4.5 5 5.5 6 

tog,o(«(H2)) [cm-3] 



6.5 



(B) 



(-70", -82") 



CO 1^ = 0.95 ± 0.05 

6-»5 



4.5 5 5.5 6 

tog,o(«(H2)) [cm-3] 



6.5 




450 - 
400 - 
350 
300 
250 
200 
150 
100 
50 



(C) 



(-60", -30") 



16.8 



CO ^ = 1.00 ± 0.05 

6-»5 



4.5 5 5.5 6 

togio(«(H2)) [cm-3] 



6.5 





450 




400 




350 




300 






250 








200 




150 




100 




50 



(D) 



(-100", 0") 



16.5 



16.6 



16.7^ 



46*^ 



CO 1^ = 0.99 ± 0.07 



4.5 5 5.5 6 

tog,o(n(H2)) [cm"'] 



6.5 



Fig. 9. Top - The gray scale and contours represent the average (logw 
scale) column density per line width (cm"- km"' s) required to repro- 
duce the observed '^CO line ratio between the peak main beam 
temperatures T,,,/, and the the peak T^i, of the '^CO 7 = 6 — » 5 line 
observed at position A (Ace = -70", AS = +32"), for different kinetic 
temperatures Tk and densities ;j(H2). Bottom - Same as top, but at posi- 
tion B {Aa = -70", A6 = -82"). 



Fig. 10. Top - The gray scale and contours represent the average (log lo 
scale) column density per line width (cm"^km"' s) required to repro- 
duce the observed '-CO ^ line ratio between the peak main beam 
temperatures T„,b and the the peak Ti„i, of the '^CO 7 = 6 — » 5 line 
observed at position C (Aa = -60", AS = -30"), for different kinetic 
temperatures Tk and densities «(H2). Bottom - Same as top, but at posi- 
tion D {Aa = -100", A(5 = 0"). 



lively. Considering a line width of ~ 6 km s we obtain sim- 
ilar column densities as in the previous two positions. That is 
A^C'^co) * 5 X 10'^ cm-2 andNi^^CO) ~3x 10'^ cm^^, for po- 
sitions C and D, respectively. At positions A and C the lines are 
optically thin, with r ranging from about 0.6 to 1 for Tk > 150 
K. At temperatures < 150 K, the lines become optically thick 
(1 < T < 3). The optically thin limit at positions B and D is 
reached at Tk ~ 120 K, with about the same ranges of optical 
depths as before, for temperatures higher or lower than 120 K. 

Assuming a density of 5x10'' cm""* and average temperatures 
of 200 K at position A and 150 K at position B, the non-LTE 
model indicates that '^^CO column densities of ~ 1.5 x 10'^ cm"^ 
and ~ 7.4 x 10'* cm"^ would be required to reproduce the ob- 
served strength (Table [1) of the '^CO 7 = 6^5 line at these 
positions, respectively. For a temperature of 200 K, the '^CO col- 
umn density at position C and D would be ~ 5.6 x 10'* cm"^. 
And the excitation temperatures would be ~ 180 K at position 
A, C and D, and ~ 140 K at position B, which are higher tem- 
peratures than estimated with the LTE approximation. However, 
these excitation temperatures are just between 10 and 20 K lower 



than the assumed kinetic temperatures, which indicates that this 
lines are close (within 10%) to the thermal equilibrium. 

From their C'**0 7 = 2 — » 1 observations, S88 estimated 
a '^CO beam-averaged column density of ~ 2 x 10'^ cm"^, 
considering a ['^CO ]/[C'**0] abundance ratio of 500. Hence, 
the column densities found for the four selected positions in 
M17 SW suggest that the warm (Tk > 100 K) and dense 
(n(H2) > 10^ cm"^) gas traced by the mid-7 '^CO lines rep- 
resent < 2% of the bulk of the cold (T^ < 50 K) and less dense 
(«(H2) ~ 10^ cm"^) gas traced by the low-7 '^CO lines. 

4.3.2. Volume filling factors 

The clump volume filling factor <t>v can be estimated from the 
ratio between the average volume density n,,,, per beam, and the 
average clump density tidump ~ 5 x 10^ cm"^ derived from the 
non-LTE model (e.g. Kramer et al. 2004). The average volume 
density per beam can be estimated from the total column density 
of the gas, and the line of sight extent of the cloud (Ddoud)- That 
is «a,, ~ A^(H2)/Dcioud- Following the work by Howe et al. (2000) 
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we can assume a '■'CO abundance ratio of 1.5 x 10 ^ relative to 
H2, and estimate the hydrogen column densities of A^a(H2) ~ 

2.3 X 10^^ cm-2, Nb(H2) ~ 1.1 x 10^^ cm^^, and A^c,d(H2) ~ 

8.4 X 10^^ cm"^, for the four selected positions. 

The line of sight extent of the cloud is a difficult parameter to 
estimate. From a 13" (~ 0. 14 pc) beam-averaged column density 
of A^(H2) ~ 8x10^"' cm"^ a volume filling factor of 0. 13 was esti- 
mated by SG90. While Howe et al. (2000) reported a ~ 0.002 
from the total column density of A^(H2) ~ 4 x 10^^ cm"^ es- 
timated at the peak column density of their '^'CO 7=1^0 
map, and assuming a cloud extent of 3 pc, which was decon- 
volved from the 4' beam of the SWAS space telescope. The line 
of sight extent should be larger than the smallest possible clump 
size (~ 0.1 pc) that we can deconvolve from our 9.4" beam. But 
we do not think it can be as large as 3 pc, which is about the size 
of the maps we present here. This holds at least for the region 
of bright '^CO and '^CO emission close to the ionization front, 
where our four selected positions are taken from. If we take the 
average between the upper (3 pc) and lower (0. 1 pc) limits of the 
cloud extent, we would obtain a cloud size of ~ 1.6 pc. This line 
of size extent of the cloud is uncertain, but perhaps more realistic 
given the geometry of M17 SW and the high resolution of our 
maps. Besides, it is similar to the diameter of the [C I] emitting 
region 1 pc) estimated by Genzel et al. (1988), and the nar- 
row spatial extension (~ 1.3 pc) of the '^CO 7 = 6^5 and [C 
I] 370 fim lines along the strip line at P.A=63° (Figures[3]&|4]l. 

Using the total column densities estimated for the four se- 
lected positions, and DcUmd = 1-6 pc, the average volume den- 
sities at position A and B would be ~ 5.3 x 10"^ cm"-' and 
~ 2.5 X 10"^ cm"-', respectively, and ~ 1.9 x 10'* cm"-' at po- 
sition C, and D. This in turn gives volume filling factors - 
nav/ndiimp of ~ 0.106, ~ 0.050 and 0.038 at positions A, B and 
C/D, respectively. These volume filling factors, as well as the 
total hydrogen densities estimated here, are larger than those es- 
timated by Howe et al. (2000), but smaller than the ones reported 
in SG90. This is an expected and reasonable result since the 
'^CO J = 6 —> 5 line traces only the warm and dense clumps and 
not the interclump medium. Besides, the volume filling factors 
estimated at the four selected positions are in close agreement to 
those estimated in other star forming regions using clumpy PDR 
models (e.g. S140, W3 Main; Spaans & van Dishoeck 1997; 
Kramer et al. 2004). 

4.3.3. Jeans stability of tine clumps 

With an average density of 5 x 10^ cm""*, and an average clump 
size of 0.2 pc in diameter, which gives a typical total clump mass 
of ~ 100 Mq in molecular hydrogen, M92 estimated that these 
clumps are not in pressure equilibrium with the interclump gas 
(with average density 3x10^ cm"^ and temperature of 200 K), 
but rather that they are self-gravitating. With these parameters, 
and a temperature of the order of 1000 K, the Jeans mass and 
radius should be about 1500 Mq and 0.3 pc, respectively. Hence 
these clumps are not near the collapsing regime. Even with our 
upper limits for the temperatures of the clumps of 230 K and 150 
K at position A and B, and 220 K and 200 K at position C and D, 
respectively, the Jeans mass and radius of these clumps would 
still be larger than those estimated with the average density of 
5x10^ cm"^. Temperatures < 170 K would be required to break 
the Jeans stability at that density. This means that the clumps at 
position B should have a slightly lower density of ~ 3 x 10^' cm"^ 
(or lower) to be Jeans stable at a temperature of about 150 K (or 
higher). 



4.4. Follow-up work 

A higher resolution map of the 609 yum (492 GHz) ^Pq 
fine-structure transition of [C I] will be obtained with FLASH on 
APEX, in order to constrain the ambient conditions of the inter- 
clump medium and the halo in M17 SW. More complex radiative 
ti-ansfer codes like RATRAN (Hogerheijde & van der Tak 2000) 
and /33D (Poelman & Spaans 2005), will be used to model the 
internal dynamics, temperature and density structure of individ- 
ual clouds. The models will also allow us to explore in detail the 
effect of absorbing foreground clouds, or multiple cloud compo- 
nents, in the line profiles. Our PDR code (Meijerink & Spaans 
2005) will provide the abundances of the molecular and atomic 
species, according to the UV flux estimated from historical data 
and our mid- J lines data. All together, these models will aid to 
test and constrain the heating and cooling of the irradiated gas. 

5. Conclusions 

We have used the dual color heterodyne receiver array of 7 pixels 
CHAMP^ on the APEX telescope, to map a region of about 2.6 
pc X 2.9 pc in the y = 6 ^ 5 and 7 = 7 ^ 6 lines of '^CO , 
the i^CO y = 6 5 and the ^Pi 370 yum (7 = 2 ^ 1) 

fine-structure transition of [CI] in M17 SW nebula. 

The completely different structure and distribution of the 
■'P2 ^Pi 370 fim emission, and its critical density, indi- 
cate that this emission arises from the interclump medium (~ 
3x10^ cm"^). On the other hand, the mid- J lines of '^CO and the 
isotope emissions, arise from the high density (~ 5 x 10^ cm"-') 
and clumpy region. 

The spatial extent of the warm gas (40-230 K) traced by the 
'^CO 7 = 7 — > 6 line is about 2.2 pc from the ridge of the 
M17 SW complex, which is smaller than the extent observed in 
the low-7 '^CO and C'^O lines reported in previous work. The 
'^CO 7 = 6^5 and [C I] 370 fim lines, have a narrower 
spatial extent of about 1.3 pc along a strip line at P.A=63°. 

An LTE approximation of the excitation temperature pro- 
vides lower limits for the kinetic temperature. The warmest gas 
is located along the ridge of the cloud, close to the ionization 
front. In this region the excitation temperature range between 
40 and 120 K. A non-LTE estimate of the ambient conditions at 
four selected positions of M17 SW indicates that the high density 
clumps (~ 5 X 10^ cm"^) cannot have temperatures higher than 
230 K. The warm (n > 100 K) and dense («(H2) > 10"* cm"^) 
gas traced at the four selected positions by the mid-7 '^CO lines 
represent ~ 2% of the bulk of the molecular gas traced by the 
low-7 '^CO lines. Volume filling factors of the warm gas rang- 
ing from 0.04 to 0. 1 1 were found at these positions. 
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